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Abstract 

We present a general framework for describing the quantum phases obtained by 
doping paramagnetic Mott insulators on the square lattice. The undoped insulators 
are efficiently characterized by the projective transformations of various fields under 
the square lattice space group (the PSG). We show that the PSG also imposes 
powerful constraints on the doped system, and on the effective action for the vortex 
and Bogoliubov quasiparticle excitations of superconducting states. This action can 
also be extended across transitions to supersolid or insulating states at nonzero 
doping. For the case of a valence bond solid (VBS) insulator, we show that the 
doped system has the same PSG as that of elementary bosons with density equal 
to the density of electron Cooper pairs. We also discuss aspects of the action for a 
d-wave superconductor obtained by doping a "staggered-flux" spin liquid state. 



1 Introduction 



It is now widely accepted that the cuprate superconductors can be under- 
stood as doped Mott insulators. Anderson pQ made the early suggestion that 
appropriate reference Mott insulating ground state should be paramagnetic 
i.e. preserve spin rotation invariance with the electron spin operators obeying 
(Sj) = on all sites, j, of the square lattice; these are also loosely referred 
to as the resonating valence bond (RVB) states of Pauling [2J, or as spin liq- 
uids. A more subtle and complete understanding of spin liquids on the square 
lattice has emerged since then. It has been found that some of the physically 
interesting spin liquids are unstable at low energies to confinement and the 
emergence of new competing order parameters. The best known examples are 
£7(1) RVB states with a gap to both spin and charge excitations, which are un- 
stable by general arguments. The instability of the simplest such state leads 
to the prominent example of competing valence bond solid (VBS) order [3], 
which will play an important role in the considerations of this paper. Other 
spin liquids can potentially be stable against symmetry breaking, [Hl5] but in- 
stead have gapless gauge and fermionic 'spinon' excitations which are strongly 
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coupled together in an 'algebraic spin liquid': we will consider the example of 
the 'staggered flux' (sF) state in this paper. The sF state also has strongly 
enhanced fluctuations of the VBS and other orders. 

Many authors have described how doping a paramagnetic Mott insulator with 
charged carriers leads to <i-wave superconductivity. [6Tr7l8ll9l[T0] This leads to a 
number of natural questions on the role of competing orders, such as VBS or 
charge density wave orders, in the doped system. (?) Can the competing orders 
have long-range order in the doped system, and how does the order parameter 
at finite doping relate to the VBS order in the undoped insulator ? (m) Does the 
long-range order co-exist with superconductivity, or can it appear only in finite 
doping insulating states ? {in) What are the theories of the quantum critical 
points between such states, and can they be of the L andau-f orbidden [TTf 1 2 II 1 3] 
variety ? (iv) In the moderate doping <i-wave superconductor, where this is no 
competing long-range order, what are the correlations of the competing order, 
and what is their interplay with vortices and quasiparticles ? This paper will 
develop a formalism designed to address such questions. We will do this in as 
general a setting as possible, relying mainly on symmetry arguments. We will 
not use any specific microscopic model of the doped antiferromagnet, and so 
will not make any quantitative predictions of the phase diagram. 

Our analysis will focus on two distinct classes of paramagnetic, undoped Mott 
insulators. In both classes, the Mott insulator is described by the dynamics 
of a compact U(l) gauge field A^ (/i = x, y,T is a spacetime index). The 
"photon" excitation of this gauge field describes S = excitations above a 
spin singlet RVB ground state. Over at least a significant intermediate energy 
range, this gauge field can be described by a conventional, Maxwell action of 
the form 

Sa = ~K cos (e^x^uAjx) ■ (1) 

□ 

We have discretized spacetime onto the sites of a cubic lattice, j, with each 
equal time section mapping onto the sites of the underlying physical square 
lattice (it should be clear from the context whether j refers to the square 
or cubic lattice), □ represents a sum over elementary plaquettes of the cubic 
lattice, A M is a lattice derivative, the indices /i, z/, A extend over x, y, r, 
and K is a coupling constant. In addition to the pure gauge fluctuations in 
Eq. ([I]), the theory of the undoped Mott insulator must also account for the 
'matter' fields: the matter associated with the average density of exactly one 
electron per site. For our purposes, it turns out that the most important effect 
[T5f5] of the matter in the undoped insulator is on its consequences for the 
monopoles in A^ (the monopoles are discrete tunnelling events in which the 
net U(l) gauge flux is changed by 2ti). The matter causes each monopole to 
transform non-trivially under the space group symmetries of the square lattice, 
and we dub these transformations (following Wen, [H] in a slight abuse of 
mathematical terminology) as the PSG (for projective symmetry group) of 
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the monopoles. [T2lfT6lfT7] The monopole PSG specifies the competing order 
parameters which emerge from the RVB state at low energies. 

We turn next to the doped antiferromagnet, where the ground state can be 
a superconductor. Here, one class of excitations of the superconductor are 
the vortices in Cooper pair condensate, each carrying electromagnetic flux 
of hc/(2e). It was shown in recent work [T6l[TT] that a general description of 
competing orders near to a superconducting state is obtained by analyzing 
the PSG of the vortices. The PSG ties each vortex to fluctuations in com- 
peting orders associated with generalized 'density- wave' and 'vorticity-wave' 
modulations, and a number of observable consequences of this connection were 
discussed. 

So how do we describe the evolution of the competing orders in the undoped 
Mott insulator to those in the superconductor ? We will argue in this paper 
that a very powerful and general method of doing this is provided by the 
PSG. Quite simply, we combine the PSG of the monopoles of the undoped 
Mott insulator with the PSG of the vortices in the doped system, and ob- 
tain a combined action which is invariant under both PSGs. In principle, this 
effective action applies in both the undoped insulator and the highly doped 
superconductors, and at all doping concentrations in between. It is designed 
to address the questions posed in the second paragraph. 

From a different perspective, the approach of this paper is an expansion about 
a "Mott quantum critical point", describing a system on the verge of super- 
conductivity due to charge condensation in a non-trivial paramagnetic Mott 
insulator. The crucial assumption is the validity of this expansion, i.e the close- 
ness in phase space of the physical system to the quantum critical point. The 
degree of stability of the non-trivial Mott insulator - the "reference state" - 
is a secondary issue. Indeed, in the first case we consider, the non-trivial Mott 
insulator is the U(l) RVB state with fully gapped spin and charge excitations, 
which is generically unstable to VBS order at low energies. We therefore de- 
note this RVB state as a VBS insulator, though most of our considerations 
never require the actual occurence of long-range VBS order of any particular 
type. Nevertheless, a number of recent works have shown that such a state 
can be a good starting point for describing the quantum critical regime in 
which VBS order is weak. The different possible paramagnetic Mott states do, 
however, generally lead to physically different predictions. A proper choice 
amongst such states must be made on energetic or phenomenological (empir- 
ical) grounds. 

The following two subsections briefly introduce the two classes of paramagnetic 
Mott insulators we consider. For the VBS insulator, introduced in Section [TTT| 
we will be able to carry our derivation of the action to completion, using the 
PSG of the monopoles described in Ref. [T2]. The sF state will be considered in 
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Section [L2l the monopole PSG of this state has not been yet been computed 
and so our analysis of this case will remain incomplete. Nevertheless, we will 
be able to make substantial progress in describing the structure of the action 
in this paper, and note significant differences from the VBS case. The final 
subsection of this introduction will briefly outline the strategy for extending 
our analysis to finite doping, and discuss the issue of flux quantization. 



1.1 The Valence Bond Solid Insulator 



Our first quantum paramagnet is obtained in a theory of quantum fluctuations 
about an antiferromagnetically ordered Neel state. This theory focusses on the 
Neel order parameter as the primary dynamic degree of freedom, and obtains 
VBS order in the state where quantum fluctuations have "disordered" the Neel 
order. 1 3 . 181 



The Neel state is described by the Neel order parameter, iij, which is a unit 
vector in spin space, related to the electron spin operator by 

Sj = Srjjiij, (2) 

where S = 1/2 is the spin, and rjj = ±1 takes opposite signs on the two 
sublattices. The action for quantum spin fluctuations of the n has a familiar 
description in terms of the 0(3) non-linear a model with Berry phases: 

5 n = 7T (A^n) 2 + iS Vj A D(n j ) ■ A T n i , (3) 

where A_o(n) = Ad(— n) is the Dirac monopole function which is used to 
yield a Berry phase proportional to the area enclosed by the worldline of each 
spin. The coupling constant g can be tune to 'disorder' the Neel state: for 
small g, the ground state has Neel order with (n) 7^ 0, while for large g, there 
are strong fluctuations of n leading to a RVB state with (n) = 0. However, 
this "disordered" state nevertheless has non-trivial quantum dynamics which 
is expressed most easily in terms of the gauge field A^ which is defined here 
by 

1 

2" 

As has been argued elsewhere, [18] for large g, the effective action for A^ in 
this RVB state is 1S4 + Sb, where 1S4 was in Eq. (TTT) and the Berry phase term 



A^ = -A D (n) ■ A M n. (4) 



S B = i2S'£A j r (5) 

3 

descends directly from the second term in Eq. ([3]). The properties of 1S4 + Sb 
have been described in detail elsewhere. [3"fT5] and will be reviewed here in 
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Section 12.11 Briefly, the Berry phase endow the monopoles in A^ with non- 
trivial transformations under the PSG, so that they transform like a VBS 
order parameter. j31T2"] Further, at low energies, the monopoles "condense" 
leading to VBS order in the ground state. The full details of the PSG will be 
presented in Section [2J but here we note the transformation of A^ under the 
time-reversal operation, T. This can be deduced from Eqs. (J3J) and (j4]): we 
have n — ► — n, and hence 

T . A T ► A T , A x * A x ) Ay ► Ay (6) 

We will observe shortly that An has a distinct transformation under T for our 
second quantum paramagnet: this distinctions plays an important role in the 
two theories of vortex dynamics in the d-wave superconductor. 

1.2 The Staggered Flux spin liquid 

This RVB state, which we denote by the sF spin liquid, has seen much recent 
discussion in the literature. |4"|5|9|2"0] and its basic properties will be reviewed 
in Section [3J Here, we note that the low energy (fixed point) properties of the 
sF spin liquid are described by Sa + S^, where describes the coupling of 
Afj, to 4 two-component, fermionic, massless Dirac fields, schematically, the 
action is S-q, = J d 2 rdrdi,[A f j], where 

U [AJ = -«*7/i + ( 7 ) 

where 7 M are the Dirac matrices. The PSG properties of the sF spin liquid are 
described in Section I3.ll Under T, the gauge field A^ transforms, as expected 
from the relativistic structure of Sq, as a conventional Lorentz vector 

T~ . A T A T , Ax Ax 5 Ay Ay. (8) 

Note the distinction of this from Eq. (jBJ) for the VBS state. Monopoles also play 
an important role in dynamics of the sF spin liquid. [2H22II14||4"] Computing 
their PSG remains an important open problem, and so we will perform our 
analysis to the extent possible, pending the eventual determination of the 
monopole PSG. 

1.3 Doping the antiferromagnet 

A detailed description of the theory obtained by doping the two paramagnetic 
Mott insulators described above appears in Sections[2]and[3]respectively. Here, 
we note some important common features of the two theories. 
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In both cases, we represent the charge carriers by two species [9|23II17] of 
spinless bosons[I] ('holons'), b\ and If the total hole density is 5, then each 
species of boson has density 5/2, and we will assume throughout that 

2 q 

where p and q are relatively prime integers. As we will see in the body of the 
paper, the number-theoretic properties of the integer q will play a crucial role 
in our analysis. 

Both b\ and b 2 carry physical electromagnetic charge e. However, they carry 
opposite charges ±1 under the internal compact U(l) gauge field A^. The 
superconducting state is obtained when the charge 2e, gauge neutral, combi- 
nation &1&2 condenses. However, with charge e bosons present in the theory, 
one may wonder if there are regions of stability of flux hc/e vortices. As we 
review below, quite generally all stable vortices have flux hc/(2e). 

Far enough from the vortex cores, we may focus on the phases of the boson 
fields alone. Let us therefore write b\ ~ e~ %dl and 62 ~ e~ td2 . We denote a 
vortex in which the phase of Q\ winds by 2n by ip\, and a vortex in which the 
phase of 9 2 winds by 2n by ip 2 . A complete dual theory of the ^1,2 vortices 
appears in subsequent sections. Here, we continue with an analysis in the direct 
picture. Far from the vortex core, the free energy contains a contribution 



(10) 



We have not included the electromagnetic magnetic field here, as its influence 
is important only at distances of order the London penetration depth. For a ipi 
vortex, V#i = cg/r, and V0 2 = at large r (eg is a unit vector orthogonal to 
the radial direction). The slow 1/r decay implies a logarithmically divergent 
vortex energy in Eq. (flU!) which is cut off at the London penetration depth. 
Let us assume that A = aeg/r at large r (which is pure gauge). Then the co- 
efficient of the logarithmically divergent term is proportional to (1 — a) 2 + a 2 . 
Minimizing w.r.t. to a, we obtain a = 1/2. Thus a ipi vortex carries a total A^ 
gauge flux = § Afj,dr M = tt. Similarly, a ip2 vortex carries A^ gauge flux of — ir. 
In both cases, at distances of order the London penetration depth, there must 
be a residual gauge flux 7r (in units of hc/e) in the physical electromagnetic 
gauge field to ultimately render the energy of a vortex finite. So both the ipi 
and ip2 vortices carry flux hc/(2e). 



1 A complementary theory with spinless fermionic holons, and bosonic spinons, is 
discussed in Ref. [24j . Such a theory is natural when the undoped VBS state is 
proximate to a deconfined critical point |11|12|13] to a magnetically ordered Neel 
state. The approach in the present paper to doping the VBS state with bosonic 
holons is more closely connected to the large ./V limit studied in Ref. [10). 
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The above A^ flux assignments also allow us to identify the connection between 
the vortices and the monopoles. A monopole is a tunneling event in which the 
A^ flux changes by 2n. Such a change in flux also occurs when a ip\ vortex 
transforms into a ip2 vortex. At low energies we expect the two tunneling 
events to always co-incide: this is especially so in the superconducting state, 
where the A^ gauge flux is confined by the condensate. So the operator 
for a monopole tunneling event is V'lV'i- 



The following Section [2] will present our theory for doping the VBS state, and 
the analogous theory for the sF state appears in Section [31 



2 Doping the Valence Bond Solid 



An explicit microscopic theory of the doping of the VBS state described in 
Section 11.11 has already appeared in Ref. [17]. There we realized the VBS 
state by a quantum dimer model, and then carried out a detailed duality 
transformation on the doped quantum dimer model. Here, we will show how 
the same results can be obtained by more abstract, but also more general, 
symmetry arguments based upon the PSG. The present analysis makes it clear 
that the results of Ref. [T7J are applicable to a far more general class of models, 
and also allows us to obtain new results on the nature of superfluid-insulator 
transitions. 



2.1 Undoped insulator 



The theory of the insulator was outlined in Section [T7TI and has been discussed 
in much detail elsewhere. Here, we focus exclusively on the crucial PSG prop- 
erties. First, we complete the mappings in Eq. (J6]) to the full set of symmetry 
operations for the underlying antiferromagnet. We will consider the following 
operations here, and in the remainder of the paper 
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T x : Translation along the x axis by one lattice site. 
T y : Translation along the y axis by one lattice site. 
R^2 '■ Rotation by 90° about a dual lattice site; 
with the origin on a dual lattice site, 
x — > y and y — > —x. 
I x nal : Reflection about the y axis of the dual lattice; 
with the origin on a dual lattice site, 
x — > —x and y — > y. 
T : Time reversal. (11) 

The mappings of are easily determined from Eqs. ([3]) and PJ, and we 
obtain 



A 

Ay 

Ay. (12) 

We have not explicitly written out the transformations of the spacetime co- 
ordinates of the fields above, because they are evident in Eq. ( FlTT) . 

The PSG of the monopoles in A^ requires careful consideration of the Berry 
phases in Eq. ([5]). As shown by Haldane, [15] each monopole acquires a net 
phase factor, which (for S = 1/2) then leads to the following PSG transfor- 
mations [12J for the monopole annihilation operator m 



T -A 

±y . 

pdual . a 
IX n/2 ■ 

rdual . A 
l x . J\ T 

T:A T 



-A • A 

-A • A 
A ■ A — 



-> -A ■ A 

-> A ■ A - 
-A ■ A - 



T x :m — > im* 



Ty-.m — > — irrv 
Jf al : m -> m 

T :m — > m. (13) 

As has been discussed elsewhere, [TTp2] e J7r//4 m has the same transformation 
properties as the VBS order parameter. The condensation of m in the compact 
U(l) gauge theory for the VBS phase then implies the appearance of VBS order 
in the ground state. 



8 



2.2 Doped antiferromagnet 



We will consider the charged excitations of the doped antiferromagnet in this 
subsection: the spin excitations will be considered later in Section 12.51 

As noted in Section 11.31 we represent the charge carriers in the doped an- 
tiferromagnet by S = 0, charge e ho Ion degrees of freedom b\ ~ e~ ldl and 
b 2 ~ e - * 6 * 2 . The boson b± carries charge +1 under A^, while the boson b 2 car- 
ries charge —1 under A^. The total density of holes is S, and the density of 
each species of boson is 5/2. 

We can now obtain an effective theory for these bosons, constrained by their 
transformations under the PSG. Interestingly, there is little arbitrariness in 
this PSG: it is almost entirely determined by the PSG of A^ in Eq. ( [121) . the 
A^ charge assignments of the bosons above, and the requirements of gauge 
invariance. We will actually not need the explicit form of the boson action; 
just the PSG will suffice for our subsequent duality mapping to the vortices. 
A simple analysis shows that the PSG of the bosons is 





■ h- 


+ & 2 ; 


b 2 - 


+ h 


T 


■h- 


+ & 2 ; 


62- 




r>dual 
K n/2 < 


:&!- 


-> ib 2 


; h- 


-> il 


rdual 

X 


■ h- 




b 2 - 




T: 


■ h- 




b 2 - 





The factors of i in the transformations under -R^ 1 are no ^ determined by 
the requirements of gauge invariance. Instead, as discussed in some detail in 
Ref . [17] , these phase factors depend upon microscopic details, and the nature 
of the short-range pairing in the superconducting ground states. The phase 
factors displayed above are those appropriate to d-wave pairing. In any case, 
such phase factors do not play any role in the duality to the vortex degrees of 
freedom. 



We proceed to apply the boson-vortex duality separately to both species of 
bosons &i 5 2- Then boson b\ dualizes to a vortex ipi ~ e - *^ 1 and a non-compact 
U(l) gauge field while b 2 dualizes to a vortex ip 2 ~ e ~ l 4>2 anc [ a non _ 
compact U(l) gauge field A 2 ^. Again, it is not necessary to write down the 
explicit form of the dual action as the PSG will suffice in determining the 
low energy continuum limit below. We can obtain the PSG of the dual fields 
by the requirements of gauge invariance, the facts that the boson currents 
A^#i — Afj. and A^9 2 + A^ must transform like the dual fluxes e^ u \ A^Aia and 
e^ lu \A u A 2 \ respectively, and conversely the dual vortex currents — 
and — A^0 2 + A 2fl must transform like the gauge flux e^xA^Ax. From these 
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requirements, it is not difficult to determine the PSG of the dual vortex matter 
fields t/)i t 2 



T 

± X 

T 
± y 

r>dual 

K -n/2 

rdual 

x 

T 



■4h 



^2 
-02 

4 
fa 



1p2 

1p2 
1p2 



4>i 

4 



(15) 



and of the dual U(l) gauge fields Ai j2 



T x : Aip 

pdual . A 
U 7t/2 ■ 71 It 



A 2m ; A 2 p 
A 2/1 ; A 2 p 



A 



2t 



Aix 





a 2t - 


-> A^ ; 


A 2a; - 


rdual 

X 


Ai T ~ 


-> -A 2r 


; 




a 2t - 


-> -A^ 


; A 2x 


T 


Al T ~ 


-> -Air 


; Ax,. 




a 2t - 


"> ~^2r 


; A 2x 



A 
A 



in 



A 2y ; Aiy 
A\y ; A 2j/ 



A 2 * 

A^ 

Ai s 
A 2 x 



.4 

.4 
.4 



l.V 

2:v 
l.v 
2;/ 



-A 2 x 
-A^ 

-> -A 



-A 



2y , 

ly 



A 
A 



ly i 

2y 



(16) 



The PSG transformations in Eqs. (fT3j) . (fT5l) . and (fl6l) contain almost all the 
information needed to obtain the vortex/monopole theory of the doped an- 
tiferromagnet. However, there is an crucial ingredient that has been missing 
from our discussion so far: the influence of the background density of 5/2 of 
each species of boson. In the dual vortex picture, this emerges as an average 
background flux of both A\ t2 of 27r(5/2) per plaquette. The consequences of 
such a background flux on the vortex PSG were studied in some detail in 
Ref. [16J, and we can simply transfer those results here separately to ipi and 
ip 2 . At a rational hole density in Eq. (jUj), the ipi j2 vortices each generate q 
degenerate low energy fields. Following Ref. [IB] , we denote these fields as <pu 
and (p 2 £ where £ = 0, 1 . . . q. So at the moment there are a total of 2q vortex 
fields that constitute the degrees of freedom of our low energy theory. The 
PSG of the these fields follow from Eq. (flBT) and Ref. [TB] 
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T x : (fu — > ^2,1+1 f2i — ► <Pi£+x 

T y : ip u — >■ <puu~ ; <^2<? — >■ V^i^^"^ 
i g-i 

• <PU ~ > —p f2m^ , 

V ^ m=0 

i g-l 

\ — 

V?2<> ~ > — 2^ Vlm 1 ^ 
\1 m=0 

^ Ual = ¥>U -> </4 ! ^2£ -> </4 

T : V9i£ -> v?i£ ; y?2£ -> ^2£, (17) 

where the £, m indices are all implicitly determined modulo q, and 

to = e 2mp/q . (18) 



In the gauge field sector, the low energy theory is most easily expressed in 
terms of the uniform and staggered dual gauge fields defined by 



Ai^ — B SfM + B ail 

A 2tl = B sf ,-B atM . (19) 



We are now ready to present the promised unification of the vortex and 
monopole PSGs. Indeed, all results follow from the PSGs in Eqs. (fl3j) and 
(TPTj) . coupled with the requirements of gauge invariance. Before proceeding, 
let us briefly discuss the latter requirements. The ipu fields have charges of +1 
and +1 under B Sfl and B afl respectively, while the ip2t fields have charges of 
+ 1 and —1. The monopole changes the A^ gauge flux by 2n, and we argued 
earlier in Section [T731 this implied that a monopole transformed like ipffli'i * n 
other words, the monopole operator m) has charges of and +2 under B Sfl 
and B afl respectively. These gauge charges are summarized in Table [H 



2.3 Low energy theory 



We are now ready to begin presentation of the effective action. This has the 
structure 

S = ( d 2 xdr \C f [A„\ + £u(i) + C 2(p + C i(p + C m ] (20) 



The first term, £^[^4^], is the action for the 5 = 1/2 fermionic spinons; we 
defer consideration of this term to Section 12.51 

The second term, Aj(i) is the action for the U(l) gauge fields, A^, B Sfl , and 
B afl . Apart from the conventional Maxwell terms, there is also a Chern-Simons 



11 





V Ul UV_>-/V 


V Ul UV_>-/V 


lVTorinnnlp 


Holon 


Holon 
















J s 


uuai 














oauge 


1 


1 


n 
U 






















Final 














Ljauge 


1 


1 

-1 


o 
z 






















Direct 












Vj (VBS) 


Gauge 








i 


-l 


and 














1 (sF) 



Table 1 



Assignments of charges under the various gauge field. The direct gauge field A^ 
is related to spin singlet fluctuations in the insulating RVB state, as discussed in 
Section [H The flux of the dual gauge field B SfM measures the electrical supercurrent. 
The dual gauge field B afl is the Chern-Simons dual of A^ i.e. the flux of B ajl 
measures the current associated with A^, and vice versa. 

term which ensures that the A^ flux is equal to the vortex current associated 
with Aip — A 2(U : this is demanded by the discussion in Section [L~3l which showed 
that there is ±n A^ flux associated with the ipi^ vortices. 




+ (^uxd u B aX ) 2 + -e^B au A x . (21) 

Further A^ dependence in the action appears only in the fermionic term 
£/[^4 M ]. If the fermions are gapped, they can be safely integrated out, at the 
cost of a renormalization of the coupling K. In this case, we can also integrate 
out A^ from Eq. (121]) . and conclude that the B afl gauge field is gapped and can 
be dropped from further considerations. The situation with gapless fermions 
is far more complicated, and requires a full analysis of the coupled A^ B afJi , 
and the fermionic excitations. 

The action C 2ip is the quadratic term for the vortex fields <pi i2 consistent with 
the gauge charge assignments and the PSG. After appropriate rescalings of 
the 921,2 fields, and rescaling of time and space co-ordinates, this has the form 
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+ s (\fu\ 2 + \ f2e\ 2 ) } + C 2ip . 



(22) 



Here we displayed the standard kinetic terms which identifies the 2q ipi^e fields 
as relativistic complex scalars. The "mass" s is the tuning parameter which 
accesses various phases of the doped antiferro magnet. For large positive s, the 
vortices are gapped, and the ground state is superconducting. Condensates 
of vortices or vortex-anti-vortex pairs can form a lower doping, leading to 
insulating or supersolid phases. 

There is an additional quadratic term, denoted £2^ above, invariant under 
all PSG transformations in Eq. (fi~7] which is allowed in some cases. For the 
special values q = 4n + 2 (n integer), it is easy to check that 



is allowed, where H is a real coupling constant. More physically, this term has 
its origin on a staggered potential that can act on the 6^2 holons. As derived 
explicitly in Ref. [17], or by the requirements of the PSG in Eq. (|14p . the boson 
Hamiltonian can contain the term 



Upon dualizing to the vortices, this term appears as a staggered flux acting 
on the vortices, which in the continuum limit generates Eq. (T23]) . 

The term is a quartic polynomial in the vortex fields (fi t 2- This polynomial 
is constrained by the PSG transformations in Eq. (|17j) . in a manner that has 
been discussed at some length in Ref. [TB]. These quartic terms control the 
precise structure of the density-wave order in the supersolid and insulating 
phases. 

Finally, we turn to the main term being introduced in this paper: the C m 
which couples the vortices and the monopoles. In obtaining this term, we 
simply have to search for invariants under the PSGs in Eq. (1131) and (j!7p . A 
straightforward analysis shows that the lowest order invariant, present only 
for q = 2n is 






3 



£=0 



+C.C, 



(24) 
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where A is a complex coupling constant. Note that this is the only term which 
couples together the tpi and y? 2 vortices, and so places a crucial role in deter- 
mining the vortex spectrum. 



2.4 Vortex spectrum in the superconductor 

The vortices are gapped in the superconductor, and here we are interested in 
the nature of the single vortex excited state. (Strictly speaking, a single vortex 
has a logarithmically divergent energy, and so the excitations are vortex-anti- 
vortex pairs — we imagine there is an anti-vortex far away, and examine the 
motion of the vortex). Once a <pi vortex has been created (say), it can trans- 
mute into a (p2 vortex by the monopole tunnelling process in Eq. f[2~4"|) . In the 
superconductor, we expect that the flux associated with this monopole will 
be largely confined to the vortex cores, and so the action for the monopole 
will be finite. Under these circumstances, it is legitimate to simply treat the 
field m as a constant [12] (in the absence of gapless fermions, it is also valid 
to replace m by a constant across the transition to the insulator). 

The spectrum of a single vortex is then simply determined by diagaonalizing 
the quadratic action £ 2 <p + £-m and ignoring the gauge fields. We describe 
the results of this procedure below for different values of q in the following 
subsections. 



2.4.1 q= 1,3 (mod 4) 

For q odd, C m = and £ 2</ > = 0, and so the tp u and ip 2 e remain vortex 
eigenstates. The degeneracy of vortex states for this case is therefore q' = 2q. 

2.4.2 q = (mod 4) 

Next, we consider the case where q is a multiple of 4. Now C m in Eq. f[2~4l) is 
non-zero, and it splits the 2q dimensional vortex space into 4x4 blocks in 
which (fu, ip2i, and v ? 2,£+g/2 are coupled to each other. In each block, 

the action has the structure 



a p\ 



tt 2 + k 2 + s + 



p a 



(25) 



a* P* 



\P* a* J 
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q 


?' 


2n + 1 


2g 


4n 




An + 2 


?/2 



Table 2 



Degeneracy of the vortex states, q', at hole density 6. Here <5 is related to q as in 
Eq. ([9|), and n is a positive integer. 

where a = A(— l) f m^e~" r,/4 , /3 = Am^e l7r//4 , f2 is the frequency, and k is the 
wavevector. The eigenvalues of this action are ± (\a\ 2 + \(3\ 2 ± (a* (3 + /?*«)) + 
Q 2 + k 2 + s. Now note that a* (3 + = 0, and hence there are 2 pairs of 
doubly degenerate eigenvalues. Consequently, the 2q vortex states have split 
into 2 blocks of q degenerate states. The lowest energy vortex states therefore 
have the degeneracy q' = q. 



2.4.3 q = 2 (mod 4) 

Finally, consider the remaining case q = 4n + 2, where n is a positive integer. 
Now C 2 ifi m Eq. (T23]) is also non-zero, but the action still has a 4 x 4 block 
structure. The matrix in Eq. ( |25|) is replaced by 



Q 2 + k 2 + s + 



( 





i^f 


a 






— i'-f 





(3 


—a 




a* 


(3* 





—i'-f 


V 


(3* 


-a* 








where 7 = (— l) e E. Now the 4 eigenvalues are pi (\a\ 2 + \(3\ 2 — ip 2 (a*(3 — (3*a)) 1 ^ 2 + 
Q 2 + k 2 + s + P27, where p\, 2 = ±1. All 4 eigenvalues are distinct, and so the 
2q vortex eigenstates have now split into blocks of q/2 degenerate states. The 
lowest energy vortex states therefore have the degeneracy q' = q/2. 



2-4-4 General discussion 

We summarize the above results on the degeneracy of the vortex states in 
Table [2j The numerology of the vortex degeneracy seems rather mysterious, 
but actually has a simple physical interpretation. The results in Table [2] can 
be reproduced by the simple formula 

^ = \- P -, (27) 
q 2 q 
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where p' and q' are relatively prime integers (as are p and q). Recalling that p/q 
is half the density of holes (Eq. (Q) away from the half- filled Mott insulator, we 
now see that p'/q' is half the density of electrons i.e. p'/q f is the total density 
of Cooper pairs per unit cell. So the degeneracy of the vortex spectrum is 
always equal to that of a model of elementary bosons with the boson density 
equal to the density of Cooper pairs. This is one of the main results of this 
paper. 

The above result can also be understood directly from the PSG. Consider the 
action of the PSG on field combinations which are neutral under B ail (only 
such combinations appear in the action above). For this, we define 

1 /2 

Tp u = ip u m 1/2 ; Tp 2l = ip 2t (rrJ) (28) 
The action of the PSG on these fields is 



T x : tpi^ — > e i7r//4 (/? 2 ^+i ; V% —>■ e m 'Vi^+i 

rp — — in/ 4 — —I — i-7i74 — 



1 



[2 ■H J U— r —pZ_^ f2m UJ 1 

V " m=0 

V* - \ E ^~ lm (29) 

V " m=0 

Note that the transformations in Eq. (129]) obey 

T x T y = —toT y T x = exp T y T x , (30) 

thus explaining the g'-fold degeneracy of vortex states. 



-tm 



2.5 Fermionic excitations 



The analysis of Ref. [T7] of spin S = 1/2, neutral excitations of the quantum 
dimer and related models applies directly to our theory of the doped VBS 
phase. These excitations are represented by Fermi operators fj a on the sites, 
j, of the square lattice, with a =f, | a spin index. This fermion carries 
gauge charge rjj, as indicated in Table [TJ Also crucial are its transformation 
properties under the PSG 
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^aa'fa', (31) 

where e aa t is the antisymmetric tensor. The transformations of /J are the Her- 
mitean conjugates of those of f a , except for the case of time- reversal. Time- 
reversal does not correspond to a canonical unitary transformation, and we 
have chosen to describe time-reversal as a symmetry of the Grassman coherent- 
state path integral. The minus sign in the time-reversal transformation of f> 
above is allowed since the p and / fields are not complex conjugates but actu- 
ally independent in the Grassman integral. The factor i in the transformation 
under Rfjf * s related to the corresponding factor in the PSG for the holons 
in Eq. (jl4| . It corresponds to an assumption of ci-wave pairing. 

With the gauge charges and the PSG at hand, we can write down allowed 
terms in the fermionic contribution, £^-[^4^], to the action: 



T x 


fa ~ 


+ u 


T 
± y 


fa- 


-fa 


7->dual 
K n/2 


fa' 


- ifa 


rdual 

X 


fa- 


-fa 


r 


fa- 


- taa'f a i 



CAM = E fl (§- - iVjAh) ho + v E 



+ J2A a e- i ^°fle cra ,fl a>(Tl +R.c. 

ja 

-t E b jl b U»,2fjafj+a,a + H.C. 



(32) 



J,a=x,y 



Clearly, A a is a pairing amplitude, and we have = — A y . 

In principle, it is now possible to take the continuum limit of Cf and then 
examine the properties of Eq. (|2Q|) . However, with the expected proliferation 
of monopoles in the insulator, we expect that fluctuations are very strong, 
and it may be more appropriate to use fields that are neutral under A^. The 
quadratic terms in the fermion dispersion in Eq. (|32l) have a gapped fermion 
spectrum, and so no low energy fermions are available to suppress monopole 
events. 

Rather, it seems more reasonable to assume that strong A^ fluctuations bind 
the bi t 2 holons to the f a fermions. We can then express the fermionic Hamil- 
tonian in terms of charge e physical electron operators Cj a : 



'■jn <'"■'■ fj„ for rjj = 1 
'> '"''•'./>, for rjj = - 



(33) 
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At the same time, we can work with the Tp 12 fields defined in Eq. ( 12 8 p . and 
then all degrees of freedom are neutral under the and B afl gauge fields: the 
action is expressed entirely in terms of the Cj a , Tp 1 2 , and B s ^ fields. 

The Hamiltonian for the Cj a fermions is essentially identical to the BCS Hamil- 
tonian 

(34) 

t ^ c^ a Cj +ayC7 + H.c. 

j,a=x,y 

Notice that the exponential in the second line containing a gauge-invariant 
combination of fields which carries physical electrical charge 2. So the c a 
fermions couple only to the physical BCS order parameter. Further, just as 
in the BCS theory, the c a fermions display the usual gapless nodal fermion 
excitations of a <i-wave superconductor. 

At this point, the coupling of the c a fermions to the vortices can be de- 
scribed following the methods discussed in a variety of papers in the literature 
[2"5]l26f27 .28 29 30J: it yields a theory for neutral fermions valid in a regime 
where the fermionic excitations can be gapless. After a singular gauge trans- 
formation which eliminates the phase of the pairing amplitude [25|26j . the 
vortices have two important effects of the quasiparticle motion: the fermions 
acquire a 'statistical' phase of n upon encircling a Tp x 2 , and also acquire a 
"Doppler shift" proportional to the local superflow velocity. The statistical 
phase is implemented by coupling both the vortices and the fermions to U(l) 
gauge fields, a M , and respectively, along with a mutual Chern-Simons term 



£ cs [cw] = - e ^v\a^d u a x . (35) 

These gauge fields have to couple to conserved currents, and for the fermions 
a convenient choice [25 .26 3 lj is the z component of the spin. We denote the 
resulting fermionic Lagrangian by £ c [a^], and refer the reader to these earlier 
works for the explicit form. For our purposes, we need the transformations 
of the a M gauge field under the square lattice symmetry operations. Using 
the usual transformations of the electron operator c a under the square lattice 
symmetry operations, and with requirements of gauge invariance, we can easily 
deduce the following PSG for 



0jl+0j + a,2+.Ajc«) / ot c 



C 'ja e aa'Cj +a , a > + H-C. 
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a x 
a y 

a y . (36) 

For the vortex sector, we couple a M to the vortex current. We assume that 
a diagonalization of the vortex spectrum has been carried out as discussed 
in Section 12.41 and focus exclusively on the q' degenerate low energy vortex 
modes. We choose these modes as eigenvectors of the T y operator, and denote 
these modes simply as Tp m , with the index m = 0, 1, . . . (q'— 1). Collecting these 
terms, our theory for the vicinity of the superconductor-insulator transition is 
then 



£dsc-VBS,v — (^fiu\d u B s \Y + C cs [a^, a M ] 

+ Ys\\( d »- iB *H ~ Kma^TPm? + S l^m! 2 f + ^ + £ c [a M ]. (37) 
m=0 ^ ' 

The action for the B Sfl field has been written in a schematic relativistic form, 
which is appropriate for short-range interactions between the bosons — the 
Coulomb interactions lead to modifications presented in Ref. [3D]. Note that 
the Afj,, .B aM , and m fields have dropped out. The term now includes 
quartic invariants in the Tp m which are invariant under Eq. (I29I) . 

The Tp m fields have charges ( m under the gauge field a M . If a M was coupling 
to the total vortex current, then we should choose all ( m = 1. However, this 
gauge field implements only a statistical phase factor of -1, and this can be 
obtained by choosing arbitrary odd integer ( m . 

An action closely related to £asc-VBS,^ has been examined previously by Lan- 
nert et al. [32] for the case of Dirac spectrum in £ c [a M ]. They considered a 
simplified model at half-filling with q' = 2. As discussed by Lannert et al, the 
Chern-Simons term in £dsc-VBS,^ drives confinement in the insulating phase 
with (Tp) ^ 0: each \I/ fermion has a a M flux tube attached to it, and this 
acquires an additional B s ^ flux tube when the total flux coupled to the Tp is 
expelled; the latter implies attachment of charge e to the spinons, and a likely 
gapping of the fermion spectrum in the insulator. A notable feature of this 
theory is that the nodal fermions remain gapless all the way up to the criti- 
cal point. A full PSG analysis on the direct lattice seems necessary to verify 
such a scenario (as in Section III.B of Ref. [16]), but below we present some 
plausible constraints under which such a critical point may obtain. 

We close this subsection by discussing some issues related to the PSG proper- 



T r : a, 



Ty . 

pdual . 
n Tr/2 ■ a r 

rdual . 
1 x ■ a r 

T : a T 



Qj t , Qj x ► & x , Qy 
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ties of £dsc-VBS,^- While the its properties should be invariant under arbitrary- 
choices for Cm-, our approximate analysis in the following subsection will find 
dependence on the values of ( m . It is therefore useful to find a choice of ( m in 
which the invariance under the PSG is as explicit as possible, without using 
non-perturbative properties of the Chern-Simons term. We will see in the fol- 
lowing subsection that it is useful to satisfy the constraint J2 m Cm = 0. This is 
clearly not possible for q' odd (which corresponds to q = An + 2), and so we 
will not consider this case further. For other values of q, we divide the vortex 
fields into 2 sets, one with ( m = 1 and the other with ( m = —1. Invariance 
under the PSG is then possible provided the set of fields with ( m = 1 either 
transform only among themselves, or all transform into fields with ( m = —1 
(and conversely for the fields with ( m = —1). For q odd (with q' = 2q), there 
is no monopole-induced mixing between the two sets of Tp u and the Tp 2l fields, 
and we can simply assign ( m = 1 for the first set, and ( m = — 1 for the second 
set. For q a multiple of 4 (with q' = q), we consider explicitly the permu- 
tative PSG for q — 4. This is specified in Eq. (C5) of Ref. [16J, and we see 
that, in this permutative vortex basis, the needed conditions are satisfied for 
Co — C2 = 1 and £1 = (3 = — 1. We suspect a similar choice is possible for other 
q multiples of 4. In all the cases for which suitable ( m are possible, the PSG 
of the Tp m implies a corresponding PSG for the a M , which we will not write 
explicitly. Combining this with the PSG for the a M in Eq. (1361) . we finally have 
to test the invariance of the Chern-Simons term /^[a^, a M ] in Eq. (|35|) : we 
find that CcJpi^a^ changes sign under most PSG transformations. However, 
this sign is clearly not physically significant because statistical phases of ±tc 
are equivalent. 



2.6 1 Undualizing ' to fractionally charged bosons 

We now examine the critical theory in Eq. ( 1371) . associated with the conden- 
sation of the Tp vortices, by undoing the duality into direct lattice degrees of 
freedom. Here it is possible to apply recent ideas [TT]fT2"] on 'deconfined crit- 
icality' to 'undualize' £dsc-VBS,^ anc ^ obtain a theory expressed in terms of 
fractionalized direct lattice degrees of freedom. The procedure for doing this 
was discussed at length in Ref. |I6] for boson models at arbitrary rational fill- 
ing on the square lattice. Here, we can apply exactly the same procedure to the 
q' vortices Tp m . As discussed in Section III. A of Ref. [16], such a transformation 
to direct lattice bosons is possible only if it is possible to find a 'permutative 
representation' of the PSG. In the following, we assume that it is possible to 
transform the relevant terms in so that such a permutative representation 
exists. Then, as shown in Section III. A of Ref. [16], the q' vortex fields Tp m 
and the U(l) gauge field B s ^ can be 'undualized' into q' boson fields £ m and 
q' non-compact U(l) gauge fields A™. The £ m boson fields each carry physical 
electromagnetic charge 2e/q', and so each Cooper pair has fractionalized into 
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q' elementary bosons. One contribution to the direct formulation of Eq. (157)1 is 
the theory of these fields presented in Eq. (3.5) of Ref. [16]. Combining these 
terms with the contribution of the Dirac fermions, we can tentatively propose 
the following 'undualized' formulation of £dsc-VBS,^ : 



g'-i 



dcconfincd 



£ \(d,-iA™)u\ +S\tr, 



m=0 



+ ]T K mn (e»u\d u Af) (e^dpAf) + E ^|^| 2 |^| 2 + C 

m,n m,n 



(3c 



where the couplings K mn and v mn can have the most general form consistent 
with the permutative PSG. Now s is the tuning parameter, which yields su- 
perducting ground states for negative values of s where the £ m bosons are 
condensed, and insulating states for positive s where the bosons are gapped. 
A specific duality transformation applied to £dsc-VBS,^, along the lines of 
Ref. [16J, does indeed yield precisely the action in Eq. fl38|) . with no residual 
Chern-Simons term — this absence is one of the principal advantages of the 
direct formulation. This duality transformation also shows that the q' gauge 
fields A™ and the gauge field are not all independent, but certain linear 
combinations are 'Higgsed' out and so become gapped. We assume that we 
are working at energies much lower than this energy gap, and so impose con- 
straints projecting out these linear combinations. In this manner we find the 
constraints 

E^ = ° ( 39 ) 

m=0 

^ = \ E tJ§. (40) 

The theory in Eq. (I38I) of q', charge 2e/q f , relativistic complex scalars £ m and 
the Dirac fermions coupled to the U(l) gauge field a M and the q' U(l) 
gauge fields A™ (subject to the constraints in Eqs. (1311 and (HU|) ). is then our 
final theory for the deconfined criticality between a (i-wave superconductor 
and VBS insulators at Cooper pair density p' /q'. 

An important ingredient so far left unspecified is the set of odd-integer val- 
ues of the ( m . A full understanding of their values probably requires a direct 
derivation of £deconfined from the underlying lattice model, as discussed in Sec- 
tion III.B of Ref. [IB] , without the long detour taken here into dual vortex 
variables. We defer such an analysis to future work, but will now note some 
important restrictions on the values of the ( m . 

To obtain these restrictions, let us consider the structure of vortices that 



21 



can be created from the £ m boson fields, using arguments similar to those 
presented in Section 11.31 Because these boson fields descended from a duality 
transformation of q' vortex fields Tp m each carrying physical magnetic flux 
hc/2e, we expect that the same vortex structure should also emerge from 
an analysis of vortex saddle points of Eq. fl38l) . That this is indeed the case 
was explained in Ref. [IB] , and we now review the argument. An elementary 
vortex is created by inducing a 2tc winding in the phase of £o (say), while 
keeping the remaining (q' — 1) boson fields £ m ^o topologically trivial. (Moving 
the phase winding to the other (</ — 1) fields yields a total of q' distinct 
elementary vortices, as expected.) Then, under the action in Eq. (13"51) . the A™ 
gauge fields will respond to minimize the total action of this vortex state. At 
a large distance, r, from the vortex center, we expect that these fields will 
be pure gauge, and oriented purely in the azimuthal (the angular co-ordinate 
9) direction; keeping in mind the symmetry of the this vortex state and the 
constraint in Eq. (1391) we obtain 

A° e (r -oo) = - 
r 

Anr^™) = -j^y ,™^0, (41) 

where A is a constant to be determined by minimizing the action. The 1/r 
decay above ensures that § dr^A™ is a constant on a contour far from the 
center of the vortex and measures the total A™ flux trapped near the center 
of the vortex. Inserting the above configurations of £ m and A™ into the action 
in Eq. (|38l) . we find that the total energy of a single vortex is logarithmically 
divergent, and the optimal vortex configuration will minimize the co-efficient 
of this logarithmic divergence. The co-efficient is proportional to [16] 

(l-^) 2 +E(W-l)) 2 , (42) 

m=l 

and minimizing this expression yields A = 1 — 1/q'. Upon including the phys- 
ical electromagnetic field, the logarithmic divergence is cutoff by the London 
penetration depth, and the total magnetic flux will be that required to sat- 
urate the remaining phase winding of the £ m at large r; a simple calculation 
from the result in Eq. (j42p shows that this yields the required flux of hc/2e. 

Turning to our objective of restricting the values of the ( m , let us now insert 
the results above for the values of A™ into Eq. (H0|) . and so obtain the value 
of the total flux associated with the vortex: 
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rn 



) 



(43) 



This total a M flux is observed by the \1/ fermions to be trapped at each vortex. 
Because the \I/ fermions pick up a -1 Berry phase around each vortex, and 
recalling the constraint that all the ( m have to be odd integers, we obtain the 
additional constraint that 



The choice ( m = 1 for all m does not satisfy this constraint. However, choosing 
an equal number of ( m = 1 and ( m = — 1 so that J2mCm = 0, does satisfy 
Eq. (jUJ); explicit choices of this type were discussed at the end of Section [2751 It 
is also clear that no solution is possible for odd q', indicating that a deconfmed 
critical point does not exist in this case. 



3 Doping the staggered flux spin liquid 

The staggered flux (sF) spin liquid has been considered in some detail in the 
context of the "SU(2) slave particle" description of Wen, Lee and collaborators 
[91120] . We will begin with the same theory in the undoped insulator, but use 
a vortex theory to describe the doped system. As noted earlier, aspects of our 
analysis will be uncomplete, because the PSG of the monopoles above the sF 
phase is not available. Our analysis will address a variety of transitions out of a 
finite doping superconducting state. The resulting quantum critical points are 
in a different universality class from the zero-doping quantum critical points 
which have been proposed elsewhere as the fixed points controlling the finite 
doping physics [33J. 

3.1 Undoped insulator 

The staggered flux spin liquid is described by a mean-field Hamiltonian ex- 
pressed in terms of a S — 1/2 fermionic spinor fj a : 




rn 



even integer. 



(44) 





(45) 
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For t' ^ t, these fermions see a staggered flux n) which apparently breaks 
translational symmetry. Since the sign of this flux cannot be removed by a 
U(l) gauge transformation, translational symmetry is implemented instead 
with a particle/hole transformation. The PSG transformations which leave 
H s p invariant are 



fj* ir l3' e ™'f}>,G 
fa -> if a 
fja irjj'Caa' fj' >a 

fja -> iVjfja ; fja ~> (46) 

As in Eq. fl3Tl) . all operations are canonical apart from time-reversal, and for 
this case the transformation of /J is independent of / CT . In all cases, the site j' 
is the image of the site j under the noted transformation e.g. j' = j + x under 
T 

Moving beyond mean field theory, the f a fermions are coupled to a U(l) gauge 
field A^, under which they have a uniform charge of +1 (see Table [1]), in 
contrast to the charge of rjj obtained in the VBS case. The PSG transformation 
of An are easily deduced from those of the f a by imposing the requirement of 
gauge invariance. The completion of the PSG Eq. (jSj), as in Eq. (fl2"|) . is 
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(47) 



Finally, to complete the analysis of the sF spin liquid, we need the transforma- 
tions of the monopoles in A^, the analog of the relations in Eq. (fT5|) . These we 
do not present here, but our analysis below can easily be extended to include 
them. 



3.2 Doped staggered flux state 



As stated earlier, we dope the spin liquid state by introducing two species of 
holon bosons, 61,2, which carry charges ±1 under the A^ gauge field. The PSG 
of these bosons can be deduced by the requirement that the physical electron 
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operator be invariant under all transformations. The latter is connected to the 
slave particles by [91TT7] 

C ja = b ]lfja + ^jb^acr'fja', ( 4 8) 

From Eqs. (H6l) and fj48|) we can deduce 
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Note that the differences from the corresponding Eq. (fT4l) for the VBS case are 
restricted to I& 2 and T. In Eq. ffl^l) . the 1, 2 vortex flavors are interchanged 
under R^jf but not under T, while in Eq. (1491) above the opposite is true. 

Proceeding with the duality from bosons b\^ 2 to vortices ipi j2 and U(l) gauge 
fields Ai^ and v42 M as in Section 12.21 the PSG of the vortices in Eq. (f!5l) now 
becomes 
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(50) 



while the PSG of the dual U(l) gauge fields A\ 2 in Eq. (fTBj) is replaced here 
by 
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Again, note that the differences from Eqs. (I15II16P for the VBS case are re- 
stricted to R^f and T. In Eqs. ( TT5lfT6l) . the 1, 2 vortex flavors are interchanged 
under R^jf but not under T, while in Eqs. (j50|5ip above the opposite is true. 
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Just as in Section 12.21 we now account for the influence of the mean hole 
density of 8 (with 8/2 = p/q) on the vortices by introducing 2q flavors of 
vortices ipu and if2i- The PSG analog of the transformations of these vortices 
in Eq. (fT7|) is now 



T 


■ <fie ~ 


-* V2,e+i i 


V>21 —> 


• <Pi,e+i 


T 
± y 


■ Vie - 




; v^f - 




?->dual 

K 7z/2 


■ <Pie - 


1 9-1 

- — E 


VlmW 






H>21 - 


- — E 


¥2mU~ 





T-.ipu^ y? 2 £ ; <^2£ — > Pie, (52) 



Once again, the differences from Eq. p7j) for the VBS case are the opposite 
treatment of the 1, 2 vortex flavors between -R^/jf anc ^ 



5.5 Lou; energy theory 



Following the analysis in Section I2.3[ we now need to write down the most 
general effective action consistent with the PSG and gauge charge assignments 
obtained in Sections 13. II and 1 3.21 This action can then describe transitions out 
of a d-wave superconductor into supersolid and insulating phases affiliated 
with the staggered flux phase. 

The required action continues to have the form in Eq. (I20I) . with the contri- 
butions Aj(i) an d C-2tp retaining their forms in Eqs. (l2Tj) and (1221) . The quartic 
term will remain unspecified as the most general quartic polynomial in ip 
which is invariant under Eq. (1521) . 

The fermionic contribution is obtained by taking the continuum limit 

of Eq. (145]) and yields the familiar Dirac form at the four nodal points, which 
was schematically indicated in Eq. (JTj). 

The additional quadratic invariant, C^, in Eq. (|22[) is now no longer given by 
Eq. ( 1231) . Instead, now there is a term whose origin is the staggered flux acting 
on the holons. Upon dualizing to vortices, this flux becomes equivalent to a 
staggered "chemical potential" acting on the vortices. A search for such terms 
reveals the following contribution which is present only when q is a multiple 
of 4: 
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(53) 



The PSG transformation associated with time-reversal allows only Eq. (1531) in 
the present case, and only Eq. (1231 in the doped VBS case. 



The final ingredient are the monopole terms C m which couple monopoles to 
vortex bilinears. As in the VBS case, the vortex bilinears must carry dual 
staggered gauge (B atl ) charge ±2, and so must mix the 1,2 vortex types. We 
will leave such terms undetermined here, but expect their influence can be 
easily included in the analysis below. 



3.4 Vortex spectrum and quantum phase transitions out of the superconductor 



As in the VBS case, it is useful to divide the discussion into various classes of 
values of g(mod 4). 



3.4.I q= 1,3 (mod 4) 

Although we do not have the explicit form of the monopole terms, C m , avail- 
able, it appears a safe assumption that such terms will not contribute for 
q = 1, 3 (mod 4). The reason for this is similar to that for the VBS case: the 
wavevectors associated with the vortex fields (fi^e are not expected to match 
those of the monopole transformations. 

The 'staggered flux' term in Eq. (|53p also does not contribute. So there are 2q 
vortex species, just as in the VBS case. The low energy theory is given by the 
sum of Eqs. (121]) . ([22"]) . (J7J) and C^. This theory will describe transitions from 
the (i-wave superconductor into proximate insulating or supersolid phases. 

The present theory is closely related to the critical theory in Eq. (j37p for 
the doped VBS case. Here the mutual statistics between the vortices and the 
Dirac fermions is implemented by B aVL gauge field, and <pi and (p% vortices 
have opposite charges under B a ^. In contrast, the mutual statistics in Eq. (1371) 
is implemented by a M , and the vortices had charges ( m . So the present case 
corresponds to choosing ( m = 1 for the (p± vortices and ( m = —1 for the 
y?2 vortices. The only remaining difference between the theories then are the 
differences between the PSGs of the vortices in Eqs. (fl7|) and (|52|) . This will 
lead to minor differences in the range of the competing orders which can 
appear in the insulating phases. 
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While the above derivation of the critical theory of fermions and vortices does 
have the advantage of preserving the lattice PSG at all stages, it does have 
the unphysical feature that the nodal points are pinned at the wavevectors 
(±7r/2, ±7r/2) even in the finite doped superconducting case. A different con- 
tinuum limit, along the lines of that discussed in Section 12.51 is needed to 
rectify this defect. 

We can undualize the critical theory to obtained the theory for the deconfined 
critical point as in Eq. fl38l) : the only changes are that the values of ( m are 
as specified above, the a M gauge field is replaced by A^ with A^ obeying the 
constraint Eq. fHUl) . and the fermion term C c [a^\ takes the Dirac from £^[^4 M ] 
in Eq. 



3.4.2 q = 2 (mod 4) 

The 'staggered flux' term in Eq. (IBUl) does not contribute for this case either. 
However, we expect non-trivial contributions from the monopole terms in C m 
now. It is a plausible hypothesis that such terms will reduce the degeneracy 
of the vortex spectrum from 2q to q/2, as was the case for the VBS state, 
reviewed in Table [2j 



3.4.3 q = (mod 4) 

This is the case most relevant for application to the cuprates, and displays 
some rather interesting features, not encountered in any of the cases considered 
so far. 

Now we have to include the 'staggered flux' term in Eq. (|53|) . This term has 
first order time derivatives, and so raises the possibility that the vortices will 
have a 'non-relativistic' dispersion spectrum. 

We also have to consider the possible influence of monopole terms, C m here. 
However, we will see below that the 'staggered flux' term in Eq. fl53|) is al- 
ready sufficient to reduce the degeneracy of the vortex species from 2q to q, as 
required by the degeneracy spectrum for this value of q from Table [2j There- 
fore, it is not an unreasonable expectation that the monopole terms will not 
contribute here. In any case, one can view the following analysis as a diago- 
nalization of the existing terms in the vortex action, which can be useful basis 
for considering the subsequent possible influence of monopole terms. 

It is useful to work with the following parameterization to diagonalize the 
vortex quadratic form in Eqs. (f22|) and (153|) . We work with the four complex 
fields W m , X m , Y m , and Z m , with m = . . . q/2 — 1, defined by 
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(54) 



Upon inserting this parameterization into Eqs. ( |22j) and ( |53j) . and ignoring 
second-order time derivative terms which are unimportant at low energies, 
the quadratic vortex Lagrangian becomes 



a/2-1 " 

E (55) 

C (J^ - iB 8T - iB aT ^j W m + -j- \(d t - iB m - i£y W m | 2 + ^ \W m \ 2 

( d \ 1 s 

f — + i£ sr + iB ar \ X m + — \{di + %B si + iB ai )X m \ 2 + j^\X m \ 2 

+ Y ™ (j^ ~ iB sr + iBar) Y m + ^ | - + ^ al )F m | 2 + ^|y m | 2 

+Z ? * n f — + i.B sr - iB aT J Z m + — \(di + iB si - iB ai )Z m \ 2 + j^\ Z m\ 2 

The full low-energy theory is now the sum of Eqs. (I2~TT) . ([7]), (|55|) and £4^. 

We note from Eq. (155!) that the fields X m , W m , Y m , and Z m are canonical 
non-relativistic Bose fields. From their B s ^ charges, we deduce that the W m 
and Y m bosons are vortices, while the X m and Z m bosons are ant i- vortices. 
There are, therefore, a total of q flavors of vortices and ant i- vortices. The 
total count of degenerate vortex/anti- vortex excitations is therefore the same 
as that obtained for the VBS case, where we also had q' = q for the case where 
q was a multiple of 4. 

However, the non-relativistic nature of Eq. (155]) has important and novel con- 
sequences for nature of the fluctuations in the superconductor. For s > 0, 
the ground state of Eq. (1551) is exactly the vacuum of the W m , X m , Y m , and 
Z m bosons. Consequently, virtual quantum fluctuations of low energy vortex- 
anti-vortex pairs are essentially totally absent in the superconductor — this is 
in strong contrast to all other cases with a 'relativistic' action, where such 
fluctations dominate and drive the superconductor-insulator transition. 

In the absence of such vacuum fluctuations, the gauge-field interactions be- 
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tween the vortices and anti-vortices are unscreened, and so will lead to the 
formation of vortex-anti-vortex bound states which will remain robust as the 
value of s is lowered. Indeed, a glance at Eq. (153]) shows that we can expect 
that the lowest energy bound states to form between the W m and the X m 
bosons and between the Y m and the Z m bosons. With the lowering of s, it is 
the energy of these bound states which will first cross zero. As the net vorticity 
of this condensing boson is zero, the this is a transition from the supercon- 
ductor to a supersolid. Thus we have reached the remarkable conclusion that 
the transition out of the sF-doped (i-wave superconductor is necessarily into 
a supersolid, for the case that q is a multiple of 4. 

We now ask whether the above spectra of low energy vortices and anti- vortices 
places any restrictions on the nature of density wave order in the supersolid. 
As in Ref. [16] we can define the density wave order operators p mn , with 
m, n = 0, 1, . . . , q — 1, at the wavevectors Q mn = (2vrp/g)(m, n) by 

3-1 

Pmn = U mn/2 [fltfU+n + <Pll<PW+n] U lm . (56) 

These operators follow from the requirement that the PSG transformations 
in Eq. (|52|) lead to transformations for p mn required for conventional density 
wave operators [16J. We now need to insert the parameterization (154)) into 
(1561) . and thence deduce the nature of p mn fluctuations under Eq. (I55|) . The 
expressions so obtained are quite lengthy, but we can understand the general 
result by considering a few representative cases. So for e.g. q = 4, we find 

p u = e -/ 4 + iW*W + X*X - iX*X x ) + (W ^Y,X -> Z) ,(57) 

while 

pio = WoAo + iW x X x + WqX* + iW{X{ + (W —>Y,X—*Z). (58) 

Notice, a crucial difference between the two cases considered above. In Eq. (1571) 
we only have combinations between W, X, Y, Z creation and annihilation 
operators; however, it is impossible to annihilate such bosons from the su- 
perconductor vacuum, and so can expect that p n fluctuations are strongly 
suppressed. In contrast, pio involves combinations of operators which are al- 
ways both creation or annihilation operators; in particular there are terms 
which lead to the creation of the W, X and Y, Z low energy bound states 
that were noted above. We therefore conclude that pio fluctuations are strongly 
enhanced as s is lowered. 

These observations can be extended into a simple general result. Let us define 
the susceptibility Xmn as the correlator of p mn and p* mn = p_ m _ n at zero 
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external frequency and momentum. Then, evaluating one-loop Xmn under the 
Lagrangian in Eq. (I55I) we find 

Xmn = 0, for m + n even, (59) 

while there is a divergent response at other values of m, n as s — ► 0: 



r dVLd 2 k 1 
Xmn 7 8vr 3 (n 2 + (k 2 + s) 2 /h 2 s ) 

~ln(A/s), for m + n odd, (60) 

with A an upper cutoff. Unless the above effects are overwhelmed by anoma- 
lously large quartic couplings in C^, we conclude that the transition from the 
doped sF <i-wave superconductor is into a supersolid in which the strongest 
density modulations are at wavevectors Q m n with m + n odd. 

We reiterate that while it is plausible that the monopoles do not modify the 
above conclusions, this has not yet been firmly established. 



4 Conclusions 



This paper has presented a new approach to the physics of doped U(l) spin 
liquids. 

We began with the two most popular examples of U(l) spin liquids on the 
square lattice at a density of one electron per site. The first, dubbed the 
VBS state, is expressed as a pure compact U(l) gauge theory; condensation 
of monopoles leads to confinement of spinons and the appearance of Valence- 
Bond-Solid (VBS) order at low energies. The second, the staggered flux (sF) 
state, has 4 species of gapless Dirac fermion spinon excitations. It is possible 
that these fermions suppress monopole condensation in the spin liquid, so that 
there is no confinement, fractionalized gapless excitations survive, and there 
is no broken lattice symmetry. Even if the sF state is unstable, it still forms 
the basis for a quantum critical point, upon which to base an effective field 
theory in the spirit of this paper. In any case, it is only more stable than the 
VBS liquid (which definitely is unstable at low energies to long-range VBS 
order) . 

Essential characteristics of the VBS and sF states are their transformation 
properties (the 'PSG') of the monopole tunnelling events under the symmetries 
of the square lattice. The existence of a non-trivial PSG for the monopole led 
to the appearance of VBS order in phases where the monopoles condense. 
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Then we doped the spin liquids with holes of rational density 5, obeying 
Eq. (Q, which defines the crucial integer- valued parameter, q. The hole degrees 
of freedom are most naturally expressed in terms of charge e, spinless bosons 
('holons'), although this does not imply existence of quasiparticles with these 
quantum numbers in any phase. The bosons appear in two species, each with 
density 5/2, carrying opposite gauge charges under the U(l) gauge field of the 
spin liquid. After a duality transformation, the charge degrees of freedom were 
encapsulated in the dynamics of vortices, each carrying magnetic flux he/ (2e). 
Central to our analysis were the PSG transformations of these vortices, and 
the couplings between monopoles and vortices that were allowed by the PSG. 

By an analysis of the low energy theories allowed by the PSG of the monopoles, 
vortices, and fermionic spinons, we arrived two main classes of results, which 
are summarized in the subsections below. 



4-1 Vortex spectrum in the superconductor 

A recurring approach in condensed matter physics is the formulation of effec- 
tive models in terms of "elementary" excitations. In some cases (e.g. Fermi 
liquids), such excitations are adiabatically connected to free electrons. Other 
elementary excitations are particular to specific states of matter, e.g. collective 
Goldstone modes resulting from symmetry breaking, or topological excitations 
such as domain walls in one dimension. The set of elementary excitations is 
usually considered to be a fundamental characteristic of a particular phase of 
matter. Most theories of quantum critical points begin by assuming the gap 
of some set of these excitations is tuned parametrically to zero, and the fields 
of the critical theory are in correspondence with them. 

Whenever gap less excitations are present, there is potential for significant am- 
biguity in the identification of elementary excitations, because new "particles" 
can be built from collections of a large number of very low energy ones, to form 
a new basis. This occurs prominently in the theory of one-dimensional systems, 
where various forms of bosonization trade gapless fermionic and bosonic field 
variables. In this paper, we sought to obtain a description of the elementary 
excitations of two dimensional superconductors on the square lattice. Clearly 
this must include both electronic quasiparticles and the vortices. When the 
former are gapless, we again have potential difficulties of the type described 
above. This paper presents one approach to deal with this formidable problem. 

The approach builds upon a previous paper, [16] in which vortices of superfluid 
phases of boson models on the square lattice were classified. A central result 
was that bosons of density p'/q' (p', q' relatively prime integers) led to vor- 
tices with a g'-fold degeneracy. Unitary transformations within this vortex 
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'flavor' space encoded various space group operators. Furthermore, a vortex 
with a fixed orientation in flavor space necessarily had static modulations at 
wavevectors (27rp'/ q') (m, n) (with m, n, integers) in all spin-singlet observ- 
ables in its vicinity. These modulations can be viewed as a strong coupling 
analog of Friedel oscillations around impurities in Fermi liquids. Quantum crit- 
icality from superfluid to insulating states can be formulated in the g-vortex 
variables. 

A "direct" attack on the analogous problem in a superconductor would be 
to try to regard the Cooper pair as the boson of Ref. [IB] . This has the dif- 
ficulty that the Cooper pair field is coupled in a very strong fashion (as in 
the Bogoliubov-de Gennes equations) to the quasiparticles. The attempt to 
extricate the collective vortex excitations from the quasiparticles in this ap- 
proach is quite non-trivial. Instead, we chose to break down the problem in 
two steps. First, we reformulated the microscopic electronic model using gauge 
theoretical methods to fractionalize the electron into a bosonic charge e holon 
and a neutral spin- 1/2 spinon. A superconducting state is then obtained as a 
superfluid state of the holons. This has the advantage that the spinons, which 
become the quasiparticles in the superconducting state, are coupled indirectly 
to the holons through a gauge field. It is then possible to treat the holons 
themselves by duality techniques directly analogous to those of Ref. [T6] . 

An interesting wrinkle in this procedure is that many different choices for the 
initial fractionalization procedure are possible. Although all of these are mi- 
croscopically equivalent, saddle point and other approximations are inevitably 
required, leading different choices to more naturally describe different insulat- 
ing states. Thus dependent upon the sort of Mott insulator reached from the 
superconductor, a particular form of fractionalization may be most appropri- 
ate. In the superfluid phase obtained by condensing the holons, however, we 
believe the same superconducting phase can be reached from different frac- 
tionalized variables. Therefore it would be natural to expect the same vortex 
excitations in the superconductor to be obtained from each of these choices. 
Despite the simplicity of this expectation, its correctness is by no means trans- 
parent in the actual calculations. Nevertheless, to the extent we are able to 
check this, we do indeed find agreement between the different gauge decou- 
plings. 

For the case of a superconductor obtained by doping a VBS state, the interme- 
diate steps of our analysis were quite involved, but our final result was simple. 
The degeneracy and PSG of the low energy vortex excitations were identical 
to that in pure boson models, with a density of bosons equal to the density of 
electron pairs i.e. the integer q' is determined here in Eq. (|27|) . Consequently 
many of the results of Ref. [16] can be applied to this electronic model without 
modification. 
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The results for the doped sF state were tentative, pending determination of 
the monopole terms. However, for the cases q = 0, 1, 3 (mod 4) we already 
obtained a vortex degeneracy identical to the VBS case in Table [2J plausibly 
suggesting that monopole terms can be ignored for these cases. It seems likely 
that a proper treatment of the monopoles would give complete correspondence 
with the VBS case, but this is left unresolved here. 



4-2 Quantum phase transitions out of the superconductor 

With the vortex spectrum in hand, we were able to address quantum phase 
transitions associated with the condensation of vortices (leading to a transition 
from the superconductor directly to an insulator) or of vortex-anti- vortex pairs 
(leading to transition to a supersolid). The critical theory for such a transition 
depends upon the fate of the fermionic S — 1/2 Bogoliubov quasiparticle 
excitations of the superconductor. In weak-coupling BCS theory, a <i-wave 
superconductor has gapless, nodal excitations at four points in the Brillouin 
zone. If these survive as gapless excitations all the way up to the quantum 
critical point of interest, then their influence has to be considered in the critical 
theory, and they could change the universality class. Alternatively, it is possible 
that the fermionic excitations are gapped at the quantum critical point: in this 
case, they can be safely integrated out and can be considered irrelevant to the 
critical theory. 

We now list the various quantum critical points obtained in our analysis in the 
following subsections. We will subdivide the discussion under headers indicat- 
ing the cases they apply to. We will explicitly write down the leading quadratic 
terms in the critical theory for each case in these subsections. Higher order 
couplings are also important, and are strongly constrained by the PSG: for 
these we refer the reader back to the body of the paper. 

4-2.1 Superconductor-insulator transition, no gapless fermions 

Doped VBS, all q: As parameters are changed in a ci-wave superconductor, 
it is possible that the electron pairing becomes short ranged, and the fermionic 
excitation spectrum is fully gapped. Such a scenario also appears plausible 
from the perspective of a confining VBS insulator that is fully gapped, as it 
is moving towards a quantum phase transition into a superconductor. In this 
case, the theory for the superconductor-insulator quantum critical point was 
found to be identical to that discussed at some length in Ref. [16] for boson 
models. Note that the notation in Ref. [16] is different from ours here — the 
integer parameter q of Ref. [16j should be set equal to the parameter q' defined 
here in Eq. ( |2"T|) . which is related to the density of electronic Cooper pairs. 
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Such a quantum critical point can occur at all values of q (defined here in 
Eq. (jUJ)) for the doped VBS spin liquid. The dual vortex theory is expressed 
in terms of q' complex scalars Tp m and a U(l) gauge field B s ^. 



Ci = ^ (e, vX d v B sX ) 2 + J2{\(d, - iB stl ) v m \ 2 + sfefj 2 } (61) 

/ m=0 L J 

Here, and in all actions below, the action for the B stl field has been written in 
a schematic relativistic form, which is appropriate for short-range interactions 
between the bosons — the Coulomb interactions lead to modifications presented 
in Ref. [30]. The superconductor-insulator transition is accessed by tuning s. 



As discussed in Ref. [16], for certain values of q' which allow a "permutative" 
PSG, the critical point can be of the 'deconfined' variety [11]: in this case it 
has a direct formulations in terms of q' complex scalars, £ m , each carrying 
electrical charge 2e/q', and coupled to q* — 1 non-compact U(l) gauge fields, 
A™. The action was given in Eq. (!38|) . but the gapped fermionic terms can be 
dropped: 



g'-i r _ 2 



m=0 



1 2 



q'-l 



+ £ K mn (e^xduAf) (^dpA^) , (62) 



m,n=0 

where the gauge fields obey the constraint in Eq. (|39l) 



IM™ = 0. (63) 

m=0 



^.^.^ Superconductor-insulator transition with gapless nodal fermions 

Doped sF, odd q: The simplest case where the nodal fermions survive 
all the way to a quantum critical point to an insulator was for the doped 
sF spin liquid, with q odd (discussed in Section 13.4.11) . In this case there is 
no non-vanishing tri-linear coupling between the monopoles and the vortices, 
and the derivation of the effective action is straightforward and preserves the 
lattice PSG at all stages. However, such a derivation has the disadvantage that 
the nodal points remained pinned at (±7r/2, ±7r/2) even in the finite doping 
superconductor. An analysis along the lines of Section [231 seems necessary to 
rectify this defect. 

There are a total of 2q vortex fields ipu and ip 2 £ coupled to 2 U(l) gauge fields, 
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B SIX and B afi , and 4 Dirac fermions, coupled to a U(l) gauge field A^. The 
critical action is 



K 2 Ks 2 K-a 2 

£3 = -g" ( £ ^a^-4a) + — {^nu\d v Bs\) + — {e^vxdyBax) 

i 9-1 f 2 

71 £=0 1 

+ 1(9^ - + iB a/t ) y^| 2 + s (\<Pu\ 2 + l^| 2 )| 
-z¥ 7M (9 M + a M )* (64) 

As above, under suitable conditions requiring the existence of a permutative 
PSG, the critical theory of Eq. AMI) can be undualized into a theory of 2q 
complex scalars £ m , each carrying electromagnetic charge e/q. The undualized 
theory is 



2?-l r 2 -, 2g-l 



m=0 



m,n=0 



- ^t, (^ + iE^-^ E^r ] *> (es) 

Z m=0 Z m=q 

where the gauge fields still obeys the constraint 

E^ = °- (66) 

m=0 



Doped VBS, g 7^ 2 (mod 4): For the doped VBS case, we have already 
discussed a situation with the absence of gapless fermions at the quantum 
critical point in Section 14.2. 1[ However, in Sections 12.51 and 12.61 we presented 
plausible conditions under which the gapless nodal fermions could survive at 
the quantum critical point in this case too. With q' defined as in Eq. (1271) . it 
was required that q' be even, and the critical theory was the same as £3 or 
£4, but with q replaced by q'/2. A more detailed direct lattice study of the 
PSG of fractionalization would be useful to firmly establish this scenario. 



Doped sF, even q: For the doped sF case, q = (mod 4), was found 
not to exhibit a superconductor-insulator transition in Section I3.4.3t rather 
a superconductor-supersolid transition obtains, which will be noted in the 
following subsection. The case q = 2 (mod 4) is likely to have monopole terms, 
and so we are not able to reach any firm conclusions for this case. 
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4-2.3 Superconductor- supersolid transition 

All the cases considered in this paper allow a superconductor-supersolid tran- 
sition: this can happen if a vortex-anti-vortex pair condensate appears before 
the condensation of single vortices. This pair condensate transforms just like 
a conventional Landau density-wave or VBS order parameter, and the critical 
theory can be developed in a traditional order parameter framework. Such 
critical theories for the superconductor-supersolid transition were discussed 
at some length in Refs. [31] and [TO] . 



Doped sF, q = (mod 4): One of our surprising results was that a 
superconductor-supersolid transition was not merely optional for a particu- 
lar case, but required as the first transition out of the translationally invari- 
ant superconductor (this conclusion is based upon an assumption on neglect 
of monopoles which has not been firmly established). This case was for the 
doped sF state with q = (mod 4). We showed in Section I3.4.3l that the £ non- 
relativistic' nature of the vortex action in this case promoted the formation 
of vortex-anti-vortex bound states which would condense first. The PSG of 
the vortex theory also placed some unusual constraints on the nature of the 
ordering in the supersolid: unless some higher order couplings where anoma- 
lously large, it was found that the supersolid in the doped sF case could have 
density modulations only at wavevectors (2irp/q)(m, n) with m + n odd. 

4-3 Discussion 

To conclude, we have found a remarkable richness in the low energy spectra 
and quantum phase transitions of a "conventional" two-dimensional <i-wave 
superconductor. The vortices of the superconducting state, when considered as 
bona fide, quantum-mechanical, quasiparticle excitations, can have surprising 
variety of wavefunctions on the lattice. These wavefunctions encapsulate the 
structure of conventional density-wave or VBS orders in proximate phases. 
Combining the quantum mechanics of these vortex quasiparticles with the 
fermionic Bogoliubov quasiparticles is a problem of considerable complexity: 
this paper has described such theories in two of the simplest cases, and the 
results are summarized above. 

We conclude by briefly discussing connections to experiments. The many 
experimental observations [3"4"|35ll36ll3"? .38 39 4 0|41II42II43] of periodic modula- 
tions in the LDOS or spin excitation spectra clearly call for attention to such 
competing orders: we have shown here, following Ref. [16], how such weak 
modulations appear naturally as an inevitable, but ancillary, consequence of 
quantum fluctuations of the vortices. Recent observations [44J of the elec- 
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tronic spectrum of such modulated states show evidence for gapless nodal 
fermions: this possibly points to proximity to deconfmed critical points with 
such fermions discussed above. For the doped sF case, we also found an un- 
expected selection of wavevectors of the modulations, which was just noted 
above: for 5 = 1/8, the density modulations were dominant at wavevectors 
(7r/8)(m,n) with m + n odd. This feature of the sF case appears to be in 
conflict with existing observations [391H2] . 
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